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Abstract 

We study charge transport driven by deformations in molecular rings and 
chains. Level crossings and the associated Longuet-Higgins phase play a cen- 
tral role in this theory. In molecular rings a vanishing cycle of shears pinching 
a gap closure leads, generically, to diverging charge transport around the ring. 
We call such behavior homeopathic. In an infinite chain such a cycle leads to 
integral charge transport which is independent of the strength of deformation. 
In the Jahn- Teller model of a planar molecular ring, X p , there is a distin- 
guished cycle in the space of uniform shears which keeps the molecule in its 
manifold of ground states and pinches level crossing. The charge transport in 
this cycle gives information on the derivative of the hopping amplitudes. 
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1 Introduction 



Adiabatic transport describes the nondissipating response of systems driven by 
time dependence ||Avron 1995] , |Kohmoto 19851 , |Kunz 1993[ [Mu 1991[ |Seiler 199 1| , 
Thouless 1994j] . An attractive feature of this theory is the geometrization of linear 
response theory. Its central result is the identification of transport coefficients with 
components of the adiabatic curvature. Because the theory describes reactive-rather 
than dissipative-transport, intuition that comes from dissipative (ohmic) transport 
is of no real value. For example, it allows for charge transport in materials that are 
nominally insulators [[Niu 1994| . The classical counterpart of the theory was studied 
in [[Berry and Robbins 1993j . 

The focus of our study will be the role of points of gap closures. This is where 
one expects interesting things to happen. This is also where one can expect to 
find universal behavior which is, to a large extent, system independent. The basic 
questions we shall ask is: What distinguishes a cycle that pinches a gap closure from 
a cycle that does not? 

The answer to this question turns out to be related to a question asked by 
Longuet-Higgins |[Longuet-Higgins 1975|| : Can one tell if a cycle of Hamiltonians 
traps a point of level crossing? The answer to Longuet-Higgins question is provided 
by looking at what is now known as Berry's phase [Perry 1984| , [Simon 1983|| . In 
the special case of Hamiltonians invariant under time reversal, the phase is ±1 
and is known as the Longuet-Higgins phase. A Longuet-Higgins phase of —1 is 
associated with conic (aka generic) level crossing, and no crossing has phase 1. The 
Longuet-Higgins phase has observable implications to the vibrational and rotational 
spectra of molecules [[Mead 1992[ [Mead and Truhlar 1979| , |Delacretaz et al. 1986 



|Whetten et al. 1986| |. What we add to this is that this phase also has observable 
implication for electronic transport. When the cycle of deformation does not trap a 
gap closure (so Longuet-Higgins phase is 1), the charge transport is proportional to 
the area of the cycle. The response to weak driving is small. In contrast, a cycle that 
traps a generic level crossing, with Longuet-Higgins phase —1, is characterized by 
homeopathic response, where the charge transport in a molecular ring is inversely 
proportional to the length of the cycle. In particular, a vanishing cause leads to a 
diverging response. In the case of infinite chains the behavior is somewhat different 
and in some ways more intriguing. It is still true that there is a weak response for 
a tiny cause when the cycle does not trap a gap closure. But, for a cycle that traps 
a gap closure, the charge transport is integral, independent of the strength of the 
driving. 

Charge pumps are quantum systems where quantized charge transport is driven 
by a periodic modulation Nju 199D| , [I'houless 1983|| , and have been recently real- 
ized [(Shilton et al. 1996| . The Thouless-Niu pumps are non- homeopathic: once the 
driving becomes weak the charge transport drops discontinuously to zero. 

Charge transport driven by deformations may be thought of as a generalized 
version of piezoelectricity. Piezoelectricity was put in a geometric framework in 
Rcsta 199^ , [King-Smith and Vanderbilt 1993j| . In this case, one replaces the global 
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question — what is the charge transport by a full cycle of deformation — by the 
local question: What is the ratio of current to the rate of deformation. Diverging 
piezoelectricity in Hubbard models of perovskites near gap closures was found in 
numerical studies in RResta and Sorella 1995fl . As we shall explain, near a generic 



crossing the ratio of current to rate of deformation is inversely proportional to the 
distance from the crossing in an infinite chain and inversely proportional to the 
square of the distance in a molecular ring. 

The divergence of transport near gap closures has a geometric interpretation. It 
is an analog of an elementary fact about the curvature of surfaces such as spheres 
and tori. The average (Gaussian) curvature (Q) of a surface with g handles, whose 
area is e 2 is, by Gauss-Bonnet, 

<«> = *lfifi>. 

The homeopathic behavior of transport is basically the divergence of the (average) 
curvature as the shrinks to zero. £ is a measure of the distance from the 

point of level crossing. In molecular rings the ratio of current to rate of deformation 
is analogous to the Gaussian curvature and diverges like s~ 2 . The charge transport 
in a cycle pinching a gap closure behaves like a line integral of the curvature and is 
of the order of 1/e. 

In the case of infinite chains with infinitely many electrons there is an extra in- 
tegration over the occupied states in the Brillouin zone. As a consequence, the ratio 
of current to rate of deformation behaves like a line integral of the curvature and 
diverges like 1/e. The charge transport in a cycle of deformation is then like an area 
integral of the curvature. For a g handle surface this gives 2(1 — g), independent 
of its area. The corresponding statement for charge transport is that a cycle of 
deformations transports an integer charge independent of the size of the deforma- 
tions (provided it encircles the point of gap closure). A model for a charge pump 
that transports integral charge for infinitesimal cycles is the insulating infinite helix 
described in section 

Quantum charge transport driven by deformation has a classical analog in the 
theory of elasticity. SX(2,R), the group of uniform shears and rotations of R 2 , is 
generated by 

^o=(_° 1 J), 9i= (I 92= (I J)- (1) 

gi and g 2 generate shears, and g generates rotations. Since 

[91,92] = 2# , 

the commutators of infinitesimal shears generate rigid rotation. Taking an elastic 
body through a small cycle of shears results in an overall rigid rotation. This 
is, perhaps, the most elementary mechanical analog which mimics the ability of 
cats to land on their feet by rotating while in free fall fShapere and Wilczek 1989b| . 
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Purcell 1977|| . Using shears to rotate is inefficient: a cycle of deformation of size e 
gives a rotation by an angle of order e 2 . The charge transported in such a cycle is 
also of order e 2 if the cycle does not pinch gap closure. If it does, the charge is of 
order 1/e. There appears to be no classical analog of this. 

A hallmark of the kind of adiabatic response discussed here is that a molecular 
quantum system can act as an ac to dc converter. For example, if one drives a cyclic 
deformation of the molecule by, say, circularly polarized radiation, the response will 
be a dc current that would give rise to a dc magnetic field. It was established 



Cohen-Tannoudji et al. 1969|| that nuclei, in his case 3 He, can also respond with a 



dc magnetic field to polarized radiation; in this case this was achieved by biasing 
the population of up and down spins. The piezoelectric response of molecular rings 
leads to the same effect by a very different mechanism. 



2 An Operator Identity 



In this section we derive the basic identity of adiabatic transport which relates the 
transport coefficients with the components of the adiabatic curvature. This identity 
is not new [[Avron 1995| , Pcllissard, van Elst and Schulz-Baldes 1994| , |Kohmoto 1985| . 
Kunz 199"3| , |JNiu 199l| , |Seiler 199~T| , |l'houless 1994| | . The purpose of this section is to 
present it as a consequence of an apparently new operator identity. 

Let H = H(x,(f)) be a family of self adjoint Hamiltonians that depend on two 
parameters x and <fi. x will stand for a deformation of the molecule and for a 
magnetic flux that threads a molecular ring. 

P = P(x, 0) is a spectral projection for H, i.e. PH = HP, with P 2 = P. 
Suppose that P is finite dimensional and smooth. This is the case, in particular, 
if H(x, <p) is a finite dimensional matrix, as in the examples we consider. The 
parameter x is slowly evolving in time in a compact interval, i.e. x(t/r) such that 
x(s) = for s < and s > 1, see fig. [Tl r is the time scale and the adiabatic limit 
is r — > oo. We shall call s = t/r scaled-time, and denote by dot derivatives with 
respect to s. The quantum evolution is described by a unitary U = U(s, <f>) solving 




Figure 1: An adiabatic switching of deformation 
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the initial value problem 



ihU = TH(x{s),<f))u, £7(0,0) = 1. (2) 

Let P = Pq{4>) denote a spectral projection for H{x(0),(j)^. p = p(s, <fi) = UP U* 
describes the physical evolution of this projection, p is always a projection, but 
it is a spectral projection only if H is time independent. Since we shall consider 
Hamiltonians that depend on time, p is not a spectral projection in general. 
The central identity of adiabatic transport says that 

rTr(pd^H) = rd^Tr(PH) + hTr(n <f)X (P)) x + O (±) , (3) 

where 

n 4tx (P) = -iP[d^p,d x P\p, (4) 

is the (fix component of the adiabatic curvature. In the time independent case, x = 0, 
this equation reduces to the Feynman-Hellmann identity. In the time dependent case 
Eq.(p]) says that the expectation value of the observable d^H at time t (in a state 
that solves the time dependent Schrodinger equation and starts as an eigenstate) is 
made of two parts: A persistent response d ( j ) Tr{PH) that survives in the er — > oo 
limit, and a linear response which is proportional to the rate of deformation x. The 
persistent response is leading in the adiabatic expansion (it is of order zero) while 
the linear response is of order 1/r. In the applications we consider d^H is the current 
operator and a; is a rate of deformation. The adiabatic curvature Q^ X (P) gives, in 
this case, the linear dependence of the current on the rate of deformation. 

This identity follows from an operator identity. For the sake of simplicity we 
formulate it under stronger conditions than necessary: 

Proposition 2.1 Let P be a spectral projection with smooth parametrical depen- 
dence on x and cfi, for a bounded self-adjoint H , so that P is separated by a gap of 
size larger than e from the rest of the spectrum. Let U be the solution of the initial 
value problem, Eq. (||) ; and p(s,(fi) = U(s,(fi)P ((fi)U*(s,(fi) with P the spectral pro- 
jection at the initial time. Then, the following operator identity holds (for all times 
s): 

r p (d^H) p = r P{d^H)P + h ^(P) x + O (—) . (5) 



The point of this identity, Eq (|5]), is that the left hand side requires the solution of an 
initial value problem while the right hand side is determined by spectral information 
alone. 

Following [ Kato 1958| and [ Avron, Seiler and Yaffe 1987 1, one introduces an aux- 
iliary evolution, XJ a = Ua(s, <fi), which respects the spectral splitting with no error. 
Respecting spectral splitting means that UaPq = PU a- A generator of such an 
evolution is: 

tH a = tH + ih[P,P]. (6) 



5 



V a solves the initial value problem: 

ihU A = TH A (x(s), ( j ) )u A , f/ A (O,0) = l. 

The evolution U A completely suppresses tunneling between the spectral subspace P 
and its orthogonal complement. The adiabatic theorem says that U A stays close to 
U: 

U A = U + 0(l/er), Wd^U = U* A d <t> U A + 0{l/er). (7) 

The adiabatic limit is re — > oo. A sufficient condition for the adiabatic theorem to 
hold is that P is smooth and the spectral gap, e, associated with P does not close. 
For precise statements and stronger results see [[Klein and Seiler 1989| , [Kato 1958 



Nenciu 1993 



From the equation of motion and the chain rule we get: 

tU* (d^H)U = ih(u*d^u). (8) 



This identity allows us to compute transport to order r _1 with error terms (which 
we do not compute) of order r -2 while using the adiabatic theorem to lowest order, 
Eq. (0). Applying Eq. (Q) to the generator of adiabatic evolution gives: 

ihu A p Q (yxd+UA) p q u* a = rp[d^H A )p = 
TP(dtH)p + ihP[P,diP\p = rPid^p + nn^iP)^ (9) 

From Eq. (|8|), and the adiabatic theorem: 

rp{d <i> H)p = r(UP U*)(d <l> H)(UP U*) = 

iKUP (U*diU)P U* = thu A p (u* A d^u A )p u* A + o(- 



= TP(dfH)p + nn^(P)x+oQ- T 

We used the adiabatic theorem in the third step. The other steps are identities. 
This completes the proof of the operator identity. d 
Tracing Eq. (H) gives: 

rTr (pd^H) = rTr{Pd^H) + hTr(^ x (P)) x + O (±) . (10) 

Since PPP = 0, 

Tr(dP)H = TrP ± (dP)PH + TrP(dP)P ± H 

= TrPP ± (dP)H + Tr(dP)P ± PH = 0, 

where P± = 1 — P is the orthogonal projection on the complement of P. This gives 
Eq. (§. 

In the adiabatic limit the persistent response dominates except in the case where 
d l pTr(PH) = 0. In this case the response in the adiabatic limit is governed by the 
adiabatic curvature. This is the case we are interested in. 

There are two mechanisms that make the persistent response vanish: 
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1. Time Reversal: Suppose that H(x,0) is time reversal invariant and that 
(d^H^j (x, 0) is odd under time reversal. Then, an elementary argument shows 
that d^Tr(PH)(x,0) = 0. There are no persistent currents in deformed 
molecules in the absence of external magnetic field, and provided time reversal 
is not spontaneously broken. 

2. Unitary Families: In the case that the <fi dependence of the Hamiltonian comes 
from a unitary transformation, i.e. 

H(x,cf>) = U(<f>)H (x)U*(<f)), (11) 

then, clearly d^Tr(PH)(x,0) = 0. Q 

Eq. (|D generalizes to the case where x is multidimensional. In the applications we 
consider x is two dimensional and naturally represented as x G (D. Assuming no 
persistent response, the adiabatic transport is: 

Tr(pd^H) = ^Re{Tr(p^(P))x} + o(±) 

= hRe {Tr(fl^(P)) ^} + O (^) . (12) 

In a molecular ring the current observable is 




where is the flux threading the ring in Gaussian units. It follows that the charge Q 
transported in a cycle in x space, in the adiabatic limit, with H (x, 0) time reversal 
invariant, is: 

Q = -c £dtTr(pd^H) = -hcRe I Tr(p^ x {P))dx 

= -e%RejTr{n^ x {P))dx. (14) 

$o — — is the quantum of flux. Note that in atomic units where e = h = 1, c = 137, 
the unit of quantum flux is a large number. This will be important subsequently, 
as we shall see. Eq. (|l4j) will be our basic tool. Note that the factors in front of the 
adiabatic curvature must be as they are on dimensional grounds. 



3 Sheared Hiickel Models 

In the Hamiltonians we consider, H(x, 0), x denotes a deformation and magnetic 
flux. The magnetic flux is virtual and is only used to define the current operator. 
The family H (x, 0) is time reversal invariant. 

lr This mechanism is responsible for the absence of persistent currents in infinite chains irrespec- 
tive of the question whether time reversal is broken or not. 
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By Wigner von-Neumann theorem, level crossings for time reversal invariant 
Hamiltonians are of codimension two. That is, like a point in a plane or a line in 
three space. For a cycle of deformation to pinch a point of level crossings, x space 
must be at least two dimensional. The simplest case is when x space is precisely 
two dimensional. 

The space of deformation of a molecule with N nuclei is 3(N — 2) dimensional 0. 
A two dimensional subspace we shall focus on is the space of uniform shears in a 
plane. 

The generators of uniform shears in two dimensions are the 2x2 matrices g± 
and g 2 of Eq. (|l]). Uniform shears can be represented by vectors x in R 2 , so that the 
shear associated to x is represented by the linear map l+xigi+x 2 g 2 - It is convenient 
to use complex notation so x is represented by a complex number, x £ (C. Under 
the shear x a point d in the sheared plane is transformed to 

d — > d + xd. 

The correspondence between the vector and complex representation is: 
x — > x = X\ + ix 2 , d — > d = di + id 2 , g\d — > d, g 2 d — > id. 
The length of d is stretched according to: 



\d\ 2 -> \d + xd\' 



\d\ 2 + xd 2 + xd 2 + \x\ 2 \d\ . 



(15) 



Suppose, for simplicity, that the hopping amplitude hij between sites i and j of 
a planar molecule is a function (typically, an exponential) of the distance \di l ' J 



| 



between the sites: h; 



K\di 



Under a small uniform deformation the hopping 



amplitudes will transform (to linear order) according to 

h(\d\ 2 ) -> h(\d\ 2 ) + 2 h'(\d\ 2 ) Re (xd 2 ). 

In particular, the Hiickel model of such a molecule is, to linear order, of the form 

H(x) = H + x(d x H) + x(dxH), (16) 

where Hq, (d x H) and (d s H) are N x N matrices. 

As an example consider the Hiickel model associated with the sheared planar 
molecular necklace X p of p sites. The tight binding model of such a necklace, see 
Fig. |, is 



H({h}) 



( 

o 





h(\d x \ 2 ) 



h(\d 2 \ 2 ) 












h(\d p \ 2 ) \ 




KK-xl 2 ] 



(17) 





\K\d p \ 2 ) . /i(|d p _i| 2 ) J 

2 The 3N nuclear coordinates include 3 center of mass coordinates and 3 Euler angles that do 
not correspond to deformations. 
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We fix the length scale so that the unstrained molecule has p equidistant atoms on 
a circle of perimeter p. The dj bond in the unstrained necklace is a complex number 
of modulus one: 



0J J 



to- 



J'-l 



10J- 



i-1/2 



\UJ 



with c<j = exp(z#), 9 = 2u/p and j = 1 . . .p. A shear x changes the distances by 




Figure 2: Molecular necklace 



Re (£x) = -Re {uo^x) 



18) 



This gives the two parameter family of Hamiltonians H(x, x). In particular 



H = h(l) 



/0 1 
1 1 
1 












1 

1 oy 



, d- x R = -h'(l) 



( 







• ^\ 


to 




























uj- 3 


U- 1 







uj- 3 ) 



(19) 



4 Adiabatic Curvature Near Crossing 



The adiabatic curvature near conic crossing was studied in [ |Bcrry 1984|| in the 
isotropic situation and by [[Simon 1983|| in the general case. The Hiickel model 
of sheared molecular rings turns out to give an intermediate situation. In addition 
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it also leads to certain non-conic crossings. We therefore explore this situation in 
some detail. 

Any two level system can be described as 

H\/%(x, 4>) = na + + no\_ + n 3 a 3 , (20) 

where our convention for Pauli matrices is: 

0\ /0 1\ f-1 



a+ = \i o)> a - = [o o) a * V 1 

A formula for the leading behavior of the x<p component of the adiabatic curvature 
is: 



Proposition 4.1 Let H(x,<p) be a Hermitian 2x2 matrix as in Eq. $iQ), with 
x G (D and <fi G R, such that to leading order in x and <p (as x,<p 0) 

n « gx m , n 3 » f<j>, / G H, <? G <D. (21) 

Then to leading order 



TrQ x </,tt± . ' , 2ra = (n + n, -i(n-n),2n 3 ) (22) 





2 


X 


2m- 2 




ri 
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Proof: A result of Perry 1984j| is that the curvature is half the spherical angle. In 



complex notation it reads: 



Trfl = ±— -— n ■ dn x dn 

2n 3 



ndn — ndn , . 

A dn 3 . (23) 



Ai \n 
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It follows that for a czrc/e of deformation, \x\ = const, the contribution of the 
adiabatic curvature to the charge transport, in the limit x — > 0, is: 



\x\ 2m [ f = 0- 

Q ± « ±e$ m7r /M'^y- « ±e$ m7r ^ |^ 2 m | x|2m ' (24) 

1^1 | J2 03 9* T U - 



Q diverges as \x\ — > if, and only if, = 0. If <f> ^ 0, then Q reaches the maximum 
value 2m7re ( i > o/(3 3 / 2 </>) for \gx m \ = \/2f(j) and vanishes in the limit. 

5 Trimers 

The smallest molecular ring is a trimer X3. H3 and Na 3 are examples. The case of a 
trimer is special in that the space of internal coordinates is three dimensional (since 
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a triangle is uniquely determined by its sides) . Three is also the dimension of linear 
transformations of the plane: Two shears and one dilation. 

A canonical choice for the internal coordinates of trimers goes back to Jacobi 



Mead 1992, Mead and Truhlar 197S]. It is instructive to relate these with the coor- 



dinates associated with uniform shears. We start by reviewing the Jacobi coordinates 



Mead 1992 



Let a, b, c be the distances between the 2-3, 3-1 and 1-2 nuclei so that a+b+c = 0. 
The Jacobi coordinates are 



w=J^a, v = J-(b-c). (25) 



Equivalently: 



a = \-u, b=\-v — \-u, c=—\—v — \-u. (26) 
V3 ' V 2 V 6 ' V 2 V 6 y J 

Evidently, when u and v are parallel the trimer is a linear molecule. The scale and 
shape of the trimer are determined by the three coordinates: 

q = u 2 + v 2 = a 2 + b 2 + c 2 , 
qX = 2 u ■ v, qY = u 2 - v 2 . (27) 

Since |X|, \ Y\ < 1 one can, without loss, set 

X = sin#cos^, r = sin#sin^, e [0, tt/2], cp G [0, 2tt). (28) 

The shape is determined by X and Y alone. This can be seen by considering the 
lengths (squared) of the sides: 

< 2 - §( 1 -5* + l/! y )- < 29 > 

Thinking of 9 and ip as the canonical spherical angles, we see that the equilateral 
triangle sits at the north pole. The equator corresponds to degenerate triangles 
where the trimer is linear. The locus of isoceles triangles are the three meridians 

tan<^ = 0, tan<£> = ±\/3. (30) 

And finally, tp — > tc — p reverses orientation. This completes the description of the 
Jacobi coordinates. 

The area (squared) of the triangle is proportional to 

A\uAv\ 2 = q 2 (l-Y 2 -X 2 ). (31) 
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We are interested on triangles which are almost equilateral, that is, near the north 
pole. To linear order in X and Y, constant area is the same as constant q. Comparing 
Eq. fl29|) and Eq. ([15]) we see that, within this linear regime, an almost equilateral 
triangle with given coordinates X and Y is obtained by applying a shear to an 
equilateral triangle, where the complex x defined in section ^ is x = ^ + i^r- F° r 
this identification we have taken an equilateral triangle oriented so that a is parallel 
to the real axis of the complex plane; rotation of the equilateral triangle by y3 /2 in 
the complex plane amounts to multiplication of x by a constant phase e~ l{f0 . 

Consider now the molecular trimer X3 threaded by a fictitious flux, so that the 
Hiickel model f\ becomes 



H(a,b,c,<j)) = a b , (32) 




where, with abuse of notation, a, b, c are shorthand for h(a 2 ), h(b 2 ), h((?). The virtual 
flux is 

£ = exp27rifi-J. (33) 

We chose a gauge so that the flux is associated with the c-bond. 
Since the characteristic equation for if (a, b, c, (f>) is 

—E 3 + E(a 2 + b 2 + c 2 ) + 2a6ccos2vr ( — ) =0, 



eigenvalues crossing occurs at 



±(a 2 + b 2 + c 2 f 2 + 3 3 / 2 a6ccos27r (J^J = 0. 

Since the geometric mean is always less than the arithmetic mean: 

(a 2 + b 2 + c 2 ) 3 ' 2 > 3 3 / 2 a6c, 

crossing can only occur if the flux is an integer or half an integer: £ = ±1. Moreover 
since the geometric and arithmetic mean coincide only if all the elements are identi- 
cal, the locus of crossings is \a\ = \b\ = \c\ and £ = ±1. The simple eigenvalue is 2a 
(the top state if a is positive) and the corresponding eigenvector is |1) = ^= (1, 1, 1). 
(Recall that £ = 1). The degenerate eigenvalue is —a and is two fold degenerate. 

We consider now a small cycle of deformations which pinches the line of level 
crossing, see fig. |3|. Such a cycle is 

3 The model is naturally associated with three (tt) electrons molecule, e.g. H3 molecule. If the 
electrons are non interacting, and spin degeneracy is taken into account, then the ground state is 
fully occupied by the two electrons and the first excited state occupied by one electron. The charge 
transported (to leading order in x) is the same as the charge transported by a single electron in 
the ground state. The currents of the two other electrons mutually cancel. 
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Figure 3: A cycle of deformation 



a(x) = h(l) + h'(l)(ux + ux); 
b{x) = h(l)+ti(l){x + x); 
c(x) = h{l) + h\\)(uJx + ujx), 



(34) 



where uj is a cube root of unity and x 6 (D runs on a small circle in the complex 
plane surrounding the origin. The phase of x has been chosen as in section |[ 

We are interested in the amount of charge moving around the triangle when the 
molecule undergoes this cycle of deformations that encircles the line of gap closures. 
In the case of a single electron, this charge can be interpreted as the total number 
of turns around the molecule that the electron makes. 

The expansion of the Hamiltonian near zero flux and near the equilateral triangle 
takes the form : 



H(x, <f)) = H + h'{\) (xV x + x V x ) + 2tt h(l) 



where 



V, 




Vrh 



(o 





-1 




















v rh 



The degeneracy splits in first order of perturbation theory, both in 
local behavior near crossing is given by the 2x2 matrix as in Eq. 



/ 



2nh{l) 



9 = 2Qh'{l), 



m 



1. 



(35) 



and in x. The 
!]]) with 

(36) 
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For a circular orbit, \x\ = const., the total charges transported per cycle in the 
almost crossing state is given by Eq. (0). If there is no external magnetic field, so 
the flux 6 = 0, one finds: 



Q=± '{7s)[Wwy (37) 

The divergence of the charge transport near conic crossing is universal for conic 
crossing. It describes the remarkable fact that the smaller the cycle that pinches 
the degeneracy, the more charge it transports. The overall constant is inversely 
proportional to the logarithmic derivative of the hopping amplitude. 

The current is not large. Only the ratio of current to the rate of driving is 
large. As the circle is shrunk, the rate of driving must also decrease in order for the 
adiabatic theory to apply. 

Driving a system by an infinitesimal cycle that pinches a crossing may or may 
not be an easy thing to do. In molecular rings, X p , this may well be a difficult 
thing to accomplish, since the Jahn Teller instability fixes a pinching cycle, or finite 
radius, which is the ground state manifold of the molecule. This point is discussed 
in more detail in the last section. 



6 Sheared Molecular Necklaces 

The study of charge transport in necklaces of p atoms with zero threading flux 
requires analysis of the gap closures for £ = 1. We take the Hamiltonian (|T9D, with 
flux dependence given by the obvious generalization of (^) . The eigenvalues of the 
unstrained molecule are: 



E m = 2 cos(m#), m = 0, 



p-1 







2tt 
P ' 



For simplicity, we deal here only with the case that p is odd. Note that the energies 
are in decreasing order, fig |j. All energies, except m = are doubly degenerate so 
that the projection on the m-th subspace is spanned by | ± m) where 



\m) = 



1 



. ,3m 



\ 1 / 



m 



0,±1, 



p 



One finds 



y/p(d s H)\m) = -ti(l)(uj m+1 +u m+1 ) 



( uo m \ 



3m+4 



V uj~ 2 J 



y/p (d<i,H)\m) 



2mh(l) 






\u m J 
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Figure 4: Eigenvalues and the roots of unity 



d s H is a double right shift: 
(m\dxH\m + k) = - 



,— ,m+fc+l i , ,m+fc+l 
h' (l ) !H ±^ (. •* ■ - 



P 



(o/W 



+ C 3fc+4 + 



-fr'(l) 5(fc + 2) {uj m ~ l + a/"" 1 ) 



Similarly, 



(mld^Hlk) = 2mh(l) 



LU k — LU m 



ATvh(l) k-m ,k + m „. 
^-uj 2 sin 9) 



(38) 



The two dimensional subspace associated with | ±m), for m — 1, . . . , is split 
in first order of perturbation theory in the <fi variable. In contrast, the split is in 
order m of perturbation theory in the x variable. All we need to compute is the 
2x2 matrix describing the crossing, and in the case £ = 1, see fig |5], we find (up to 
an overall scale): 



H 1/2 = 2(-/i'(l))' 



/"X 1 cos(-m + 2k - 1)9 N 

\k=i cosm9 — cos(m — 2k)9 / 

4tt/i(1) . , . 
— - — smlmV) a 3 <p. 

p<S>o 



cos m 



(m-l)0 (pVfV+ + (oA) m a_) 



(39) 



This puts us in the framework of proposition |4.1| with explicit values for g, m and 
/. Note that as far as degenerate perturbation theory is concerned, there are also 
other terms that split the degeneracy, like x<p. However all these are irrelevant for 
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ill- 




Figure 5: Coupling of states by the perturbation 

computing curvature since g, m and / are all we need in Eq. fl22f) . The transport is 
given by Eq fl24|). 

We want to conclude this section with two cautionary remarks about the scope of 
validity of the results of this section. The moral of both of these is that the results 
for conic crossings have a wider scope of validity than the results for non-conic 
crossings. 

1. The theory of deformed tight binding Hamiltonians was developed to first 
order in x. For gaps that open to first order, that is for conic crossings, 
such a perturbation expansion is consistent. However, for gaps that open in 
higher orders, one should have considered sheared tight binding models valid 
to higher orders in x. To lowest order, Eq. fl39|) may or may not remain 
applicable, depending on whether there exist paths that connect \m) to | — to) 
more efficiently than in Fig. [5]. 

2. We have restricted the space of deformations to uniform shears. When other 
deformations are also possible, as is the case for planar molecules X p with 
p > 4, these cannot be neglected: it is always possible to build a deformation 
that will open the gap in first order and will therefore provide the lowest order 
perturbation. 
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7 Crossing Bands and Direct Integrals 



In this section we describe a version of the equation for adiabatic transport, Eq (|||), 
for systems described by direct integrals. This is the case for infinite periodic chains 
of non- interacting electrons that admit Bloch decomposition ||Grossmann 1972] . From 



a physical point of view there are two important consequences to this extra struc- 
ture. The first is that charge transport becomes integral. The second, and related, 
is that the formulas for charge transport, which for finite molecular rings involved 
non-universal factors (the logarithmic derivative of the hopping amplitudes), be- 
come truly universal. Studies of Chern numbers in the context of charge transport 



have been made in 


Bcllissard, Bovicr and Chez 1991 




Kunz 1993, Thouicss 1983 




Kreft and Seiler 1996, Novikov 1981 


, Poilblanc et al. 1991, Tan 1994 


, Wilkinson 1984 


Yakovenko 1991 





Consider the Hamiltonian represented by a direct integral (over the torus): 

H B (x) = 



H x, 



dcj) 
2^' 



H(x, <p) acts on a fixed Hilbert space, say 7i. In this and the following section 
stand for Bloch momentum. A spectral projection on an energy band is 



will 



PbM 



P(x,4>), the spectral projection for H(x,<ft), is smooth and periodic in <p. Consider 
the observable 

(d^H B )(x) = f dtH(x,<j>)^. 

This is interpreted as current: It is made up from velocity d ( / ) H(x, <fi) and density |2. 
In this case, <fi runs on the unit circle. The persistent response for infinitely many 
non interacting fermions filling the band vanishes by periodicity: 



Tr (PBd^Hs) = j* ^ Tr (PH) = 0. 



(40) 



Nearly Crossing Bands: The question we address is what can one say about 
"axial" integrals of the adiabatic curvature 



J < ^TrVt cj , x [P(x, 



-/, I —Tr(P(x, 
'-Tr 2tt 



[d^P(x,<P),d x P(x,<P)}), (41) 



in the case of nearly crossing bands. More precisely, suppose that x = is a band 
closure for some value of <fi. What is the analog of Eq. (|24|) in this case? 

Suppose that P(x, <fi) is a projection in C 2 associated with the two level Hamil- 



tonian Eq. (pUP, with n(x,(p) and n^(x,x,(f)) as in proposition |4.1| and we set the 
origin of <p space at the zero of n and 713, assuming approach to zero as in Eq. (pip. 
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In general, there may be several such zero points. The contributions of these be- 
have additively and it is therefore enough to consider one. Since the d<p integral is 
dominated by the singularity: 



A similar formula holds for f^(P). Now the divergence is like x~ l and the power 
is independent of the order of crossing. The order shows up as a linear factor. 

In the case that there are several such points the adiabatic response collects the 
contribution from all of them. These come with one sign for the top gap closure and 
the opposite sign for the bottom gap closure, fig. |6|. 



E 




1 1 — phi 

Pi 2 Pi 

Figure 6: Top and bottom gap closures 

For fixed Fermi energy in a gap, with one point of gap closure, the total charge 
transport by a loopf] encircling a gap closure, is precisely : 

Q = Tm. (43) 

Homeopathy means that integral charge is transported even if the deformation 
is minuscule. 



± 



on i 



n- 



4i(|n| 2 + n§) 3 / 2 dx2ir 



m\g\ 
m 



\x 



2m- 2- 



-X 



Wl\X 

m 



-oo {\gx m \ 2 + \fcj)\ 2 fl 2 
+ 0(1) 



+ 0(1) 



Anix 



-oo (1+0 2 ) 3 / 2 

+ 0(1). (42) 
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Figure 7: Shearing a helix in the transverse plane 

8 Homeopathic Charge Pump 

Consider an infinitely long helical chain: A single pitch of the helix is made of an 
open circle with p atoms per pitch. The cross section of the helix is a circle and the 
unstrained atoms on the circle are equidistributed. Of course, there is no reason in 
the world why a chain of identical atoms would want to form a helix. So we may 
imagine that this structure is enforced by an external constraint. We shall consider 
uniform shears of the helix in a plane perpendicular to the axis of the helix. 

A tight binding model of such a helix is described by a p x p matrix H{x,<f>), 
where x G C is the deformation and x = is the undeformed helix. <fi stands for 
the Bloch momentum, traditionally denoted by k. H(x, (f>) is identical to the tight 
binding (Hiickel) Hamiltonian of a necklace of p atoms enclosing a flux <p. Such 
Hamiltonians were analyzed in section ^|. 

Consider a quantum system of infinitely many (non interacting, spinless) elec- 
trons on such a helix so that there are q electrons per p atoms. The unstrained 
chain, x = 0, has continuous spectrum in the interval [—2, 2] without any gaps. It 
is convenient to think of this as a collection of p bands with closed gaps: 

[-2, 2 cos(tt(1 - 1/p))}, [2 cos(tt(1 - 1/p)), 2 cos(tt(1 - 2/p))], ... [2 cos(n/p), 2}. 
4 Not necessarily a circle. 
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Shearing the chain in the plane perpendicular to the chain axis, see fig. [7], will break 
the spectrum into p bands, so that the lowest q bands will be occupied. Consider 
now a cycle of deformations \x\ = e which avoids the origin. Such a cycle pinches 
the gap closure at x = 0. (See section |[) Collecting results from the sections below 
gives the following result: 

A small cycle of shears in the plane perpendicular to the axis of an infinite helix, 
with pitch p (p odd), transports integer (non zero) charges from — oo to oo even if 
the orbit is arbitrarily small (but not zero). The amount of transported charge is 
q/2 if q, the number of electrons per pitch of the helix, is even. If q is odd, then the 
charge transport is —(p — q)/2. 

These results can be seen as follows. By Floquet theory, the band edges are 
determined by H(x, 0) and H(x, it). The two cases are simply related when p is odd 
since the Hamiltonian, Eq. fllTD , satisfies 

H({h},0)U=-UH({h},ir), (44) 

with a fixed unitary 

Using this, the spectral and geometric results for the zero flux of section ^| can be 
easily translated to the = 7r case. 

The lowest gap near 2cos(8(p — l)/2) opens in (p — l)/2-th order of degenerate 
perturbation theory. By symmetry, Eq. (fP|), this then also holds for the top gap. 
The second gap from the bottom opens, likewise, as the second gap from the top. 
This opens in first order of degenerate perturbation theory, etc. This is explained in 
fig. f)| The general rule which gives the order in which a gap opens depends on the 
parity of the gap. The results are summarized in a table that also gives the Chern 
number of the gap: 



Gap 


Order 


Chern 


2j 


j 


j 


2j+l 







9 Elastic and Magnetic Jahn Teller Instabilities 

The analysis carried out so far assumed that the deformation x and the flux are 
externally controlled parameters which can be varied at will. In reality, both are 
dynamic variables with their own equations of motion. The question one needs to 
consider is how this may affect transport. 

For example, in our analysis it is important that time reversal is not sponta- 
neously broken and the flux was a virtual flux. The dynamical equations could, in 
principle, force ^ 0. If this was to happen then homeopathic response is censored. 

Even for a molecular trimer X 3 (in three dimensions), the full dynamics is a 
formidable problem which brings out all the intricacies of the three body problem 
Mead 1992| , |IVlead and Truhlar 1979] |Englman 1972] |Bersuker 1984] . Similarly, the 



full dynamics of is determined by QED. 
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The static problem is much simpler. It is determined by an energy functional 
whose minimizers are classical equilibrium configurations. For the nuclear coordi- 
nates the minimizer gives the ground state configurations and for the flux it gives 
Ampere's law. Such a functional is : 

£(x,0) + Q(x) + ^0 2 . (45) 

Here E(x, <fi) = E(x, <p + <3?o) is an eigenvalue of say the Hiickel Hamiltonian, Q(x) is 
a classical potential energy and L is a geometric constant (L/c 2 is the self-inductance 
of the molecule). Note that the period of E(x, <fi) in <ft is large since in atomic units 
$o — c = 137 is large. 

The minimizer of this functional with respect to the flux gives the analog of 
Ampere's equation (in Gaussian units) 

o<p L c L 

Let us now specialize to the case of planar molecular rings X p and let us assume 
that the elastic energy Q(x) depends only on the distances \dj\ between neighboring 
atoms. From Eq. flTS|) the elastic energy of a molecular necklace X p , in the Harmonic 
approximation, is an isotropic function of x: 

QW = fE»)|-K(0)|) 2 = 7 ^(u^x+u^x) 2 



2 
pK 



x\\ (47) 



where K is the "spring constant". From Eq. (|39|), the electronic energy surface of 
the one-electron state near crossing is also isotropic in x, and the energy functional 
is of the form: 

P K t ,2 , 1 



± sj\g\ 2 \x\ 2m + (/0) 2 + ^\x\ 2 + — 2 . (48) 

In this functional g, K, L are all of order one but / is a small number of order l/$o- 
It is convenient to use rescaled variables 

V = x m , <^ = 0/$ o . (49) 
In these variables the energy functional is 



P K l.A2/m , ^ 



2 



± ^\g\ 2 \y\ 2 + (f%)v 2 + F —\y\ 2/m + (so) 

Here g, K, L and /$o are ai l °f order one in atomic units, while $o — 137. 
The minimizers of this functional follows: 

1. With the positive sign for the square root the minimizer is the trivial solution 

(f) = {p = y = x = 0. 



21 



2. With the negative sign for the square root and conic crossing, m = 1, the 



minimizer has 



(p = and \x\ = ^ = . = because its "spring 

constant" is stiff, of order 10 4 . This solution is the classical Jahn Teller in- 
stability. Because the energy functional is isotropic, the minimizer is a circle. 
The charge transported in this cycle is 



Q = ± 



n 2 eKh(l)sm(2n/p) 
2 [h'{l)] 2 



(51) 



3. With the negative sign for the square root and non-conic crossing, m > 2, the 
stiff term in floDp is now the elastic spring, which is not harmonic. Therefore the 
minimizer is x = and \<p\ = L\f\. This is a magnetic Jahn Teller instability, 
where there is spontaneous breaking of time reversal. Time reversal is only 
weakly broken because / is small: the spontaneous flux is of the order of 
10 _4 $o- As we have already stressed in a previous section, conclusions drawn 
from sheared Hiickel models in situations where the crossing is not conic, are 
fragile and depend sensitively on perturbations. This may be a reason why 
magnetic Jahn Teller has not been observed. 




Figure 8: The preferred cycle of deformation. 
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